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Particle  Tracking  Model  Transport 
Process  Verification:  Diffusion  Algorithm 

by  David  B.  King  Jr.  and  Tahirih  C.  Lackey 


PURPOSE:  This  Dredging  Operations  and  Environmental  Research  Technical  Note  (DOER- 
TN)  describes  a  new  set  of  test  cases  developed  to  investigate  the  behavior  of  the  diffusion 
algorithm  contained  within  the  Particle  Tracking  Model  (PTM).  The  test  cases  were  designed  to 
examine  different  aspects  of  diffusion  in  isolation  and  in  combination.  The  test  cases  needed  to 
be  complex  enough  to  mimic  real-world  conditions  but  idealized  enough  to  have  analytical 
solutions  for  comparison.  These  test  cases  were  developed  to  achieve  two  objectives:  verification 
of  the  diffusion  algorithm  and  the  establishment  of  benchmark  test  cases  to  confirm  that  future 
model  updates  will  not  have  inadvertently  degraded  capabilities.  This  technical  note  describes 
how  the  model  outputs  compare  with  analytical  predictions. 

INTRODUCTION:  The  PTM  is  a  Lagrangian-based  particle  tracking  model  developed  by  the 
U.S.  Army  Corps  of  Engineers  (US ACE)  Engineer  Research  and  Development  Center  (ERDC)  to 
efficiently  track  large  numbers  of  particles  as  they  move  through  complex  flow  fields.  It  allows  the 
modeler  to  predict  deposited  and  suspended  sediment  densities  in  space  and  time  along  with  final 
particle  fates  (Demirbilek  et  al.  2004;  Davies  et  al.  2005;  McDonald  et  al.  2006;  Lackey  and 
McDonald  2007).  Although  a  versatile  model  currently  utilized  in  various  coastal,  estuarine,  and 
riverine  applications,  PTM  is  specifically  designed  to  predict  the  fate  of  material  suspended  during 
dredging  and  placement  operations  and  to  address  the  stability  and  fate  of  in-place  sediment 
including  dredged-material  mounds,  sediment  caps,  and  contaminated  sediment  deposits.  PTM 
models  the  physical  processes  of  advection,  diffusion,  settling,  deposition,  burial,  and 
resuspension,  as  appropriate  for  each  parcel  at  each  numerical  time-step,  to  simulate  the  transport 
of  a  specified  sediment  distribution.  Thus,  required  PTM  inputs  include  a  characterization  of 
both  the  bed  sediments  and  the  tracked  sediments.  Instead  of  undertaking  the  impossible  task  of 
modeling  every  grain  of  sand,  silt,  and  clay,  sediment  is  discretized  into  parcels.  Each  parcel  is 
representative  of  a  specific  mass  of  sediment.  In  aggregate,  the  specified  sediment  parcels 
preserve  the  overall  size  distribution  and  total  mass  of  the  sediment  source.  Along  with  position 
through  time,  each  parcel  is  tagged  with  a  variety  of  additional  attributes  such  as  mass,  density, 
grain  size,  and  suspension  status. 

PTM  requires  the  input  of  hydrodynamics  (i.e.,  water  surface  elevation  and  velocities),  defined 
upon  a  bathymetry  grid  that  is  provided  through  an  external  model.  This  allows  the  modeler  to 
perform  multiple  PTM  runs  to  compare  alternative  sediment  releases  exposed  to  the  same 
hydrodynamic  conditions  without  the  computational  overhead  of  regenerating  flow  conditions 
for  each  hydrodynamic  run.  PTM  operates  within  the  Surface-water  Modeling  System  (SMS), 
which  provides  effective  visualization  tools,  making  it  useful  for  assessment  of  dredging 
practices  and  proposed  dredging  operations. 


Approved  for  public  release;  distribution  is  unlimited. 


Form  Approved 
OMB  No.  0704-0188 


Report  Documentation  Page 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

JUL  2015 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2015  to  00-00-2015 

4.  TITLE  AND  SUBTITLE 

Particle  Tracking  Model  Transport  Process  Verification:  Diffusion 
Algorithm 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Army  Engineer  Research  and  Development  Center  „Vicksburg„MS 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  Dredging  Operations  and  Environmental  Research  Technical  Note  (DOER-TN)  describes  a  new  set  of 
test  cases  developed  to  investigate  the  behavior  of  the  diffusion  algorithm  contained  within  the  Particle 
Tracking  Model  (PTM).  The  test  cases  were  designed  to  examine  different  aspects  of  diffusion  in  isolation 
and  in  combination.  The  test  cases  needed  to  be  complex  enough  to  mimic  real-world  conditions  but 
idealized  enough  to  have  analytical  solutions  for  comparison.  These  test  cases  were  developed  to  achieve 
two  objectives:  verification  of  the  diffusion  algorithm  and  the  establishment  of  benchmark  test  cases  to 
confirm  that  future  model  updates  will  not  have  inadvertently  degraded  capabilities.  This  technical  note 
describes  how  the  model  outputs  compare  with  analytical  predictions. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

13 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ERDC  TN-D0ER-D18 
July  2015 

PTM  DIFFUSION:  PTM  uses  a  random  walk  model  to  calculate  the  diffusion  velocity.  The 
horizontal  component  of  the  dispersive  velocity  (  Ud )  is  assumed  to  be  isotropic  and  is  computed  as 


Ur 


(m/s) 


(1) 


where  the  horizontal  eddy  diffusivity  (  £d)  is  defined  as 


eD  =  KEthu  (m2/s) 


(2) 


the  shear  velocity  ( u  )  as 


* 

u 


(3) 


and  the  bottom  roughness  (ks)  as 


ks  =  3  *  Dgo  (m)  (4) 

In  these  equations,  77  is  a  random  number  uniformly  distributed  between  0  and  1,  U  is  the  free 
stream  flow  velocity  (m/s),  KEt  is  the  horizontal  turbulent  diffusion  scalar,  h  is  the  flow  depth 
(m),  and  D90  is  the  90th  percentile  sediment  diameter  (i.e.,  90%  of  the  sediment  grains,  by 
weight,  have  diameters  smaller  than  this  value)  (m).  These  equations  are  all  discussed  in  further 
detail  in  MacDonald  et  al.  (2006)  and  Lackey  and  MacDonald  (2007). 

The  vertical  component  of  the  diffusion  velocity  (Wd)  is  computed  similarly  as 

ILD=2(n-0.5)^  (m/s)  (5) 

Based  upon  Fischer  et  al.  (1979),  PTM’s  vertical  eddy  diffusivity  (sy)  algorithm  has  been 
updated  from  that  given  in  MacDonald  et  al.  (2006)  to 


Sy=-KEJJ 


z(h-  z)2 


(m2/s) 


(6) 
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where  KEv  is  the  vertical  turbulent  diffusion  scalar  and  z  is  the  particle’s  vertical  position  in  the 
water  column  (m).  It  is  seen  that  the  vertical  eddy  diffusivity  has  a  parabolic  dependence  upon 
the  vertical  location,  and  at  mid-depth,  Eq.  6  reduces  to 


sv=KEvhU/ 16  (m2/s)  (7) 

PTM  allows  the  modeler  to  scale  the  turbulence  levels  in  the  model  by  accepting  user-defined 
values  for  KEt  and  KEv.  One  of  three  turbulence  conditions  (sd  =  0,  £y=  0,  or  sD  =  sv)  was  used  in 
each  of  the  tests  in  this  report.  The  third  condition  was  satisfied  by  releasing  neutrally  dense 
parcels  at  mid-depth  and  setting: 


K 


Ev 


l6g  y 

hU 


16KEtu 

U 


(8) 


METHODS:  The  series  of  test  cases  described  in  this  report  were  designed  to  examine  diffusion 
along  each  axis  independently  and  jointly.  Tests  1-7  employed  a  steady  one-dimensional  (ID) 
flow  field  that  was  aligned  with  the  axis  of  the  channel  (i.e.,  the  x-direction).  Tests  1-3  looked  at 
diffusion  in  the  x-  (down-channel),  y-  (cross-channel),  and  z-  (vertical)  directions,  respectively, 
in  isolation.  Tests  4-6  examined  pair-wise  diffusion  in  the  x-  and  v-,  the  x-  and  z-,  and  the  y-  and 
z-  directions,  respectively.  Additionally,  Test  6  looked  at  y-  and  z-directed  diffusion  using  two 
different  parcel  release  protocols  (Tests  6a  and  6b).  Test  7  examined  diffusion  occurring  along 
all  three  axes  simultaneously.  Tests  8  and  9  also  examined  diffusion  in  all  three  axis  directions. 
Test  8  used  an  unsteady  ID  flow  (a  combined  sinusoidal  and  steady  flow  that  resembled  aspects 
of  a  tidal  flow  in  a  river).  Test  9  employed  a  steady  two-dimensional  (2D)  flow  field  (rotational 
flow  about  a  central  point  that  resembled  aspects  of  a  vortex  flow).  Table  1  lists  the  conditions 
for  each  test. 

For  all  the  tests,  PTM  released  neutrally  buoyant  parcels  at  mid-depth  in  a  flat  basin.  Flows 
experienced  zero  friction  (slip  conditions)  at  all  side  and  bottom  boundaries.  That  is,  flows  did 
not  include  logarithmic  boundary  layer  profiles.  In  addition,  the  tests  were  designed  such  that 
parcels  remained  well  away  from  side  and  bottom  boundaries  at  all  times. 

Each  of  the  tests  employed  one  of  four  types  of  parcel  release  protocols:  instantaneous  cross- 
flow  line  source  release,  instantaneous  along-flow  line  source  release,  instantaneous  point  source 
release,  or  continuous  point  source  release.  Examples  of  these  four  types  of  releases  are  shown  in 
plan  view  in  Figure  1,  panels  A-D,  respectively.  Each  of  these  examples  shows  the  release  point 
as  a  black  line  or  point  centered  near  the  left-hand  (upstream)  end  of  the  panel  and  a  cloud  of  red 
parcels  at  a  time  after  release  or  after  initial  release.  The  clouds  of  parcels  have  expanded  over 
time  by  diffusion  as  they  have  been  advected  downstream  (i.e.,  in  the  x-direction)  from  the 
release  site. 
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Table  1.  Test  cases. 

Test 

Name 

Axes  of 
Diffusion 

Source 

Type 

Horizontal 

Diffusion 

Vertical 

Diffusion 

Flow  Type 

1 

X 

ICfLS 

On 

Off 

1 D  Steady 

2 

y 

IAfLS 

On 

Off 

1 D  Steady 

3 

Z 

IPS 

Off 

On 

ID  Steady 

4 

*,y 

IPS 

On 

Off 

ID  Steady 

5 

X,Z 

ICfLS 

On 

On 

ID  Steady 

6a 

y,z 

IAfLS 

On 

On 

1 D  Steady 

6b 

y,z 

CPS 

On 

On 

1 D  Steady 

7 

x,y,z 

IPS 

On 

On 

1 D  Steady 

8 

x,y,z 

IPS 

On 

On 

ID  Sinusoidal 

9 

x,y,z 

IPS 

On 

On 

2D  Steady 

IPS  =  Instantaneous  Point  Source  (Figure  1-C) 

ICfLS  =  Instantaneous  Cross-flow  Line  Source  (Figure  1-A) 
IAfLS  =  Instantaneous  Along-flow  Line  Source  (Figure  1-B) 
CPS  =  Continuous  Point  Source  (Figure  ID) 


RESULTS 

Test  1.  This  test  examined  PTM  diffusion  in  the  x-direction  in  isolation  that  is  generated  by  a 
ID  steady  flow.  The  analytical  solution  of  the  ID  diffusion  equation  (Fischer  et  al.  1979)  is 


C{x,t)  = 


M 

v/4^y 


exp 


(kg/m) 


(9) 


where  C(x,t)  is  the  ID  particle  concentration  (kg/m),  M  is  the  total  mass  of  particles  released 
(kg),  and  t  is  time  (s). 

This  concentration  can  be  converted  to  a  density  (kg/m  )  by  assuming  a  constant  distribution  of 
this  concentration  that  is  1  m  high  by  1  m  wide.  The  analytical  standard  deviation  (a)  for  this 
Gaussian  distribution  is 


cr  =  ^2sj  (m)  (10) 

The  location  of  the  centroid  (=  Ut)  is  a  function  of  time  in  the  x-direction.  Its  position  remains 
unchanged  at  its  midchannel  location  in  y  and  at  its  mid-depth  location  in  z. 

For  this  PTM  run,  parcels  were  released  instantaneously  along  a  cross-flow  line  (Figure  1A). 
This  type  of  release  eliminated  the  effects  of  diffusion  in  the  cross-flow  (y )  direction  by  making 
it  constant  so  long  as  the  parcel  sampling  scheme  stayed  well  clear  of  the  ends  of  the  line. 
Parcels  experienced  no  vertical  motion  because  gravity  was  nullified  by  creating  neutrally  dense 
particles,  and  vertical  diffusion  was  eliminated  by  setting  Kev  =  0. 
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A 


B 


# 


C 


D 


Figure  1 .  Examples  of  four  types  of  parcel  release  protocols.  See  text  for  explanation. 

The  effects  of  along-flow  directed  diffusion  were  measured  by  examining  the  distribution  of 
parcels  contained  within  the  black  (solid  line)  rectangular  sampling  region  shown  in  Figure  2A. 
The  cross-flow  limits  on  the  sampling  region  were  set  by  having  them  be  >  5  cr  from  the  centroid- 
translated  ends  of  the  release  line.  In  the  along-flow  direction,  the  sampling  region  was  centered 
about  the  analytical  centroid  of  the  distribution  (i.e.,  Ut),  and  the  width  was  set  at  ±4.59<r.  The 
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sampling  region  was  then  divided  into  5 1  bins,  each  extending  the  width  of  the  sampling  region 
in  the  cross-flow  direction  and  being  0.1  8  <7  wide  in  the  along-flow  direction.  Figure  2B  shows  a 
blowup  of  the  area  of  Figure  2A  within  the  dashed  green  box  and  includes  a  section  of  the 
centroid-translated  release  line  (blue  line)  and  the  ends  of  the  5 1  sample  bins  (black  lines). 

This  binning  procedure  was  applied  in  a  similar  way  for  all  the  tests.  By  trial-and-error,  it  was 
determined  that  this  binning  methodology  provided  an  appropriate  resolution  of  underlying 
distribution  without  overly-smoothing  the  data.  The  number  of  parcels  within  each  bin  was 
counted,  and  the  sum  was  converted  to  a  concentration.  A  comparison  of  the  PTM  parcel 
distribution  with  the  expected  Gaussian  distribution  is  shown  in  Figure  3,  panel  A. 


Figure  2.  Sampling  region  and  binning  scheme.  See  text  for  discussion. 


Figure  3.  Comparison  results  for  Tests  1  and  2. 


Comparisons  of  the  results  for  all  the  tests  are  shown  in  Table  2.  In  this  table,  the  individual  tests 
are  in  the  different  columns.  A  comparison  of  x-directed  diffusion  created  by  PTM  with  that 
predicted  by  the  analytical  solution  is  given  in  the  upper  block  of  rows,  y -directed  diffusion  in 
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the  center  block  and  z-directed  diffusion  in  the  lower  block.  “St  Dev  Ratio”  (standard  deviation) 
is  the  PTM  calculated  values  divided  by  the  analytical  values.  The  analytical  values  for 
“Skewness”  and  “Kurtosis”  are  zero.  “Max  Ht  Ratio”  is  a  ratio  of  the  maximum  heights  of  the 
two  distribution  curves  (PTM/ Analytical).  “Area  Ratio”  is  a  comparison  of  the  areas  under  the 
two  distribution  curves  (PTM/ Analytical),  and  “Corr  Coef’  (correlation  coefficient)  is  a 
statistical,  pair-wise  comparison  of  the  values  comprising  the  two  curves. 


Table  2.  Comparison  statistics  between  measured  and  analytical  results. 

|  Test  # 

1 

2 

3 

4 

5 

6a 

6b 

7 

8 

9 

#  of 

samples 

66,640 

37,261 

24,576 

47,120 

100,000 

15,843 

St  Dev 
Ratio 

1.00092 

1.05231 

1.00141 

1.01248 

0.95754 

0.97175 

Skewness 

-0.00064 

0.00787 

-0.00509 

0.00501 

0.00084 

-0.00896 

X 

Kurtosis 

-0.12560 

-0.01067 

-0.20061 

-0.06432 

0.03118 

0.13095 

Max  Ht 
Ratio 

0.99299 

1.00293 

0.98751 

1.02431 

1.10306 

1.08147 

Area 

Ratio 

0.99951 

1.00000 

0.99997 

0.98688 

1.00000 

0.99545 

Corr  Coef 

0.99957 

0.99984 

0.99958 

0.99921 

0.99780 

0.99464 

1 

#  of 

samples 

66,615 

38,290 

24,741 

3,872 

47,104 

100,000 

16,410 

St  Dev 
Ratio 

0.96771 

1.02275 

0.95071 

1.00792 

0.99841 

0.95340 

0.98120 

Skewness 

0.00154 

-0.00018 

-0.00173 

0.04556 

-0.00167 

0.00072 

0.09682 

Y 

Kurtosis 

-0.01725 

-0.23361 

0.13691 

0.15142 

-0.12758 

-0.09190 

0.06496 

Max  Ht 
Ratio 

1 .02892 

0.96449 

0.96966 

0.98702 

0.98814 

1.01129 

1.06047 

Area 

Ratio 

0.99913 

1.00000 

0.93020 

0.96295 

0.98655 

1.00000 

1.03108 

Corr  Coef 

0.99885 

0.99949 

0.99463 

0.99395 

0.99944 

0.99942 

0.99771 

1 

#  of 

samples 

100,000 

25,028 

26,634 

3,786 

45,081 

100,000 

16,363 

St  Dev 
Ratio 

0.97025 

1.05216 

1.06527 

0.99666 

0.96051 

0.91991 

0.91293 

Skewness 

0.00049 

0.00106 

-0.00657 

-0.02119 

-0.00163 

0.00061 

0.02048 

z 

Kurtosis 

-0.11943 

-0.00570 

0.03067 

-0.02953 

0.04459 

-0.11803 

-0.36331 

Max  Ht 
Ratio 

1.02159 

0.96718 

0.98798 

1.00167 

1.02985 

1.10069 

1.03160 

Area 

Ratio 

1.00000 

0.99997 

1.00138 

0.95219 

0.94418 

1.00000 

1.02812 

Corr  Coef 

0.99969 

0.99933 

0.99897 

0.99541 

0.99805 

0.99774 

0.99773 

Test  2.  This  test  examined  PTM  diffusion  in  the  y-direction  in  isolation  that  is  generated  by  a 
ID  steady  flow.  The  influence  of  diffusion  in  the  x-direction  was  eliminated  by  making  it 
constant  in  a  way  analogous  to  that  described  for  Test  1.  The  instantaneous  along- flow  line 
source  type  of  parcel  release  is  shown  in  Figure  IB.  Figure  3B  compares  the  model  output  with 
the  analytical  solution. 
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Test  3.  This  test  examined  PTM  diffusion  in  the  z-direction  in  isolation  that  is  generated  by  a 
ID  steady  flow.  The  influence  of  diffusion  in  the  x-  and  y-directions  were  negated  by  setting  KEt 
=  0.  Parcels  were  released  as  an  instantaneous  point  source.  At  a  later  time,  they  would  still 
appear  as  a  point  in  plan  view  (rather  than  as  a  cloud  as  shown  in  Figure  1C).  In  side  view,  the 
distribution  of  parcels  would  appear  as  a  vertical  line  (overlapping  row  of  dots)  centered  at  mid¬ 
depth  that  increases  in  length  as  it  is  advected  downstream.  Figure  4 A  compares  the  model 
output  with  the  analytical  solution. 


Test  4.  This  test  examined  PTM  diffusion  in  the  x-  and  v-  directions  generated  by  a  ID  steady 
flow.  Vertical  diffusion  is  turned  off.  Parcels  were  released  as  an  instantaneous  point  source  as 
shown  in  Figure  1C.  The  analytical  concentration  per  unit  area  (Fischer  et  al.  1979)  is 


C(x,y,t)  = 


M 

4m;nt 


4sDt 


>  (kg/m2) 


(11) 


The  standard  deviation  of  this  distribution  is  also  given  by  Equation  10.  Figure  4B  compares  the 
model  output  for  this  test  with  the  analytical  solution. 

Test  5.  This  test  examined  PTM  diffusion  in  the  x-  and  z-directions  generated  by  a  ID  steady 
flow.  Parcels  were  released  as  an  instantaneous  cross-flow  line  source  (Figure  1A),  and  cross- 
flow  diffusion  was  eliminated  as  described  for  Test  1.  However,  for  this  test,  the  vertical 
diffusion  was  turned  on.  Equation  1 1  describes  the  analytical  concentration  (with  the  y  variable 
replaced  by  the  z  variable).  Figure  5A  compares  the  model  output  with  the  analytical  solution. 

Test  6a.  This  test  examined  PTM  diffusion  in  the  y-  and  z-directions  using  an  along-flow 
instantaneous  line  source  release  (Figure  IB)  generated  by  a  ID  steady  flow.  This  test  is  similar 
to  Test  2  with  the  vertical  diffusion  turned  on.  Figure  5B  compares  the  model  output  with  the 
analytical  solution. 
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Test  6b.  This  test  provided  an  alternative  release  methodology  to  examine  PTM  diffusion  in  the 
y-  and  z-directions  from  that  presented  in  Test  6a.  This  test  examined  vertical  and  horizontal 
cross-flow  transects  through  a  plume  (Figure  ID).  The  steady  state  analytical  solution,  from 
Fischer  et  al.  (1979), 


C(x,y,z) 


M 


4tT£dX 


exp 


4 £dx 


(kg/m2) 


(12) 


Figure  5.  Comparison  results  for  Tests  5  and  6a. 


is  valid  provided  that  the  transects  are  sufficiently  far  upstream  from  the  head  of  the  plume.  In 
this  equation,  Ms  is  the  mass  rate  of  parcel  release.  The  standard  deviation  for  this  solution  is 


g  -  yjls^xlu  (m)  (13) 

For  this  test,  the  PTM  output  was  analyzed  at  a  position  >  5  crupstream  of  the  head  of  the  plume 
(defined  as  the  Ut  position,  with  to  being  the  time  of  first  release).  Figure  6A  compares  the  model 
output  with  the  analytical  solution.  The  increased  level  of  noise  in  the  PTM  solutions  for  this 
case  is  largely  due  to  the  smaller  numbers  of  parcels  within  the  sampled  region  (Table  2). 

Test  7.This  test  examined  PTM  diffusion  in  the  x-,  y-,  and  z-directions  generated  by  a  ID  steady 
flow  with  an  instantaneous  point  source  release  (Figure  1-C).  The  analytical  solution  (from 
Fischer  et  al.  1979)  is  given  as 


C(x,y,z,t)  = 


M 


(4  nsDt) 


3/2 


exp 


-\x-JJif  -y2  -z2 
4sDt 


(kg/m3) 


(14) 


The  standard  deviation  of  this  distribution  is  again  given  by  Equation  10.  Figure  6B  compares 
the  model  output  with  the  analytical  solution. 
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Figure  6.  Comparison  results  for  Tests  6b  and  7. 


Test  8.  This  test  examined  PTM  diffusion  in  the  x-,  y-,  and  z-directions  generated  by  a  ID 
unsteady  flow.  The  parcels  were  released  as  an  instantaneous  point  source  (Figure  1C).  The  flow 
was  primarily  sinusoidal,  but  the  addition  of  a  small  steady  (x-directed)  component  meant  that  at 
the  end  of  each  cycle,  the  centroid  was  advected  farther  downstream.  For  this  test,  the  eddy 
diffusivity  is  no  longer  a  constant  but  changes  with  each  model  time-step  because  of  the 
changing  velocity,  as  seen  in  Equations  2  and  6.  To  calculate  the  analytical  concentration 
(Equation  14)  and  standard  deviation  (Equation  10)  for  this  comparison,  an  root-mean-square 
value  of  the  velocity  was  used  to  calculate  the  eddy  diffusivity.  For  convenience,  the  time  picked 
for  a  comparison  was  at  the  end  of  a  flow  cycle.  Figure  7  A  compares  the  model  output  with  the 
analytical  solution. 


Figure  7.  Comparison  results  for  Tests  8  and  9. 


Test  9.  This  test  examined  PTM  diffusion  in  the  x-,  y-,  and  z-directions  generated  by  a  2D 
steady  flow.  The  horizontal  flow  field  rotated  about  a  central  point,  mimicking  solid-body 
rotation.  Parcels  were  released  as  an  instantaneous  point  source  at  a  position  away  from  the 
center  of  the  flow  and  analyzed  at  the  end  of  one  revolution  (Figure  8).  A  graphical  comparison 
of  the  PTM  distributions  with  the  analytical  solution  is  shown  in  Figure  7B. 
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SUMMARY  AND  CONCLUSIONS:  A  series  of  tests  were  conducted  to  verify  the  proper 
functioning  of  the  diffusion  algorithm  in  PTM.  These  tests  examined  diffusion  along  all  three 
axes  independently  and  jointly.  Most  tests  used  a  ID  steady  flow  field,  but  some  used  more 
complex  flows,  including  ID  unsteady  flow  and  2D  steady  flow.  All  tests  compared  the 
distribution  of  particles  generated  by  PTM  with  their  analytical  solutions.  The  results  of  the 
comparisons  are  provided  in  Table  2.  All  of  the  types  of  comparisons,  including  the  correlation 
coefficients,  show  strong  agreement  between  the  PTM  particle  distributions  and  the  analytical 
distributions.  Therefore,  for  all  the  test  cases  examined,  this  research  found  that  the  PTM 
diffusion  algorithm  behaves  as  expected. 


Figure  8.  2D  flow  grid  and  parcel  positions  after 
one  complete  flowfield  revolution  for 
Test  9. 
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